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ABSTRACT 

We study the emission by dust and stars in the Large and Small Magellanic Clouds, a pair of low- 
metallicity nearby galaxies, as traced by their spatially resolved spectral energy distributions (SEDs). 
This project combines Herschel Space Observatory PACS and SPIRE far- infrared photometrjF] with 
other data at infrared and optical wavelengths. We build maps of dust and stellar luminosity and mass 
of both Magellanic Clouds, and analyze the spatial distribution of dust/stellar luminosity and mass 
ratios. These ratios vary considerably throughout the galaxies, generally between the range 0.01 < 
^dust/i* < 0.6 and 10~ 4 < Ma ust /M* < 4 x 10~ 3 . We observe that the dust /stellar ratios depend 
on the interstellar medium (ISM) environment, such as the distance from currently or previously 
star- forming regions, and on the intensity of the interstellar radiation held (ISRF). In addition, we 
construct star formation rate (SFR) maps, and find that the SFR is correlated with the dust /stellar 
luminosity and dust temperature in both galaxies, demonstrating the relation between star formation, 
dust emission and heating, though these correlations exhibit substantial scatter. 

Subject headings: galaxies: general - infrared: galaxies - galaxies: ISM - dust, extinction - galaxies: 
evolution - Magellanic Clouds 



1. INTRODUCTION 

Among the nearby galaxies, the Small Magellanic 
Cloud (SMC) and Large Magellanic Cloud (LMC) rep- 
resent unique astrophysical laboratories for studies of 
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the lifecycle of the interstellar medium (ISM), because 
of their proximity (« 50 and 60 kpc for the LMC and 
SMC, respectively; Ngeow & Kanbur 2008; Szewczyk et 
al. 2009), low metallicity (w 1/5 and 1/2 Z Q ; Dufour 
et al. 1982; Russell & Dopita 1992), and favorable view- 
ing angle (24° for the LMC; Nikolaev et al. 2004). The 
SMC and LMC (MCs) are well-suited for detailed stud- 
ies of galaxy evolution processes, either as whole galax- 
ies or on the scales of individual star-forming regions or 
resolved stellar populations. In particular, these galax- 
ies are suited for studies of the ISM phases, interstellar 
dust properties, and star formation rates as a function of 
metallicity and radiative environment (i.e., the intensity 
and hardness of the radiation field). 

The MCs are also unusual as relatively bright 'satellite' 
galaxies (compared to other satellites) very close to the 
Milky Way (MW). The luminosity distributions of satel- 
lite galaxies in groups have been the subject of many sta- 
tistical studies (e.g., Skibba et al. 2007; van den Bosch et 
al. 2007), and some recent studies have focused on MW- 
like galaxies with MC-like satellites in the Sloan Digital 
Sky Survey, Galaxy And Mass Assembly, and in mod- 
els and simulations (e.g, Liu et al. 2011; Tollerud et al. 
2011; Sales et al. 2011; Busha et al. 2011; Robotham et 
al. 2012). In addition, analyses of the MCs' star forma- 
tion, structures, and velocities have yielded evidence for 
tidal interactions between the LMC and SMC: the Mag- 
ellanic Bridge and SMC Wing appear to consist of tidally 
stripped material; there is tidally triggered star forma- 
tion within the MCs and the Magellanic Bridge; and the 
LMC's offset and warped bar are arguably due to a previ- 
ous collision with the SMC (Yoshizawa & Noguchi 2003; 
Mastropietro et al. 2005; Harris 2007; D'Onghia & Lake 
2008; Gordon et al. 2009, 2011; Bekki 2011; Besla et al. 
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2012; Cignoni et al. 2012). 

Nearly half of the bolometric luminosity of the Uni- 
verse is channeled through the mid- and far-infrared (IR) 
emission of galaxies (e.g., Hauser & Dwek 2001). There- 
fore, in order to develop a more complete picture of 
galaxy evolution, it is necessary to study the IR emis- 
sion that is reprocessed by dust in galaxies. We seek to 
understand the physical processes that regulate galaxy 
evolution, including the formation of stars and the inter- 
stellar radiation field (ISRF), and the return of radiant 
energy from these stars into the ISM. 

IR observations of galaxies have yielded insight into the 
nature and spatial distribution of their dust content, in- 
cluding its extinction and emission properties (e.g., dust 
temperature), as well as its composition, grain size dis- 
tribution, and abundance (e.g., dust mass). A variety of 
studies have investigated these issues in the MCs, often 
utilizing mid- and far-IR observations, in order to con- 
strain their dust properties (e.g., Weingartner & Drainc 
2001; Bot et al. 2004, 2010; Sakon et al. 2006; Bernard 
et al. 2008; Paradis et al. 2009, 2011; Gordon et al. 2010; 
Sandstrom et al. 2010; Galliano et al. 2011; Planck Col- 
laboration 2011). By combining IR and other data, some 
studies have also explored how dust production, destruc- 
tion, and heating are related to atomic and molecular 
gas content, star formation, and the interstellar radi- 
ation field. For example, constraints on the dust and 
gas content are sufficient to estimate gas-to-dust ratios 
and deviations from its expected metallicity dependence 
(Stanimirovic et al. 2000; Gordon et al. 2003; Leroy et 
al. 2007; Galliano et al. 2011). 

In general, analyses of the spectral energy distributions 
(SEDs) of galaxies yield information about the relation 
and balance between dust and stellar emission, by ex- 
ploiting multi- wavelength observations (e.g., Draine et 
al. 2007; Dale et al. 2007; da Cunha et al. 2010; Skibba 
et al. 2011). The spatial distribution of dust vis-a-vis stel- 
lar luminosity and mass is the main focus of this paper. 
The distribution of dust/stellar ratios within galaxies, 
and its connection to star formation and dust heating, 
has important implications about the energy conversion 
from stellar to dust emission. 

It is now possible, with data from IRAS, Spitzer, and 
Herschel, to constrain the radial distributions of dust and 
stars within galaxies (Sauvage et al. 1990; Zaritsky 1999; 
Munoz-Mateos et al. 2009; Mattsson & Andersen 2012), 
and to analyze the spatial distribution of dust, gas, and 
starlight in the MW (Paradis et al. 2012), M31 (Smith 
et al. 2012), M33 (Boquien et al. 2011), and M83 (Foyle 
et al. 2012), and other nearby galaxies (Aniano et al. 
2012; Galametz et al. 2012). Others have added optical, 
Ha, and UV observations to constrain the evolution and 
environmental dependence of star formation in the LMC 
and SMC (Bell et al. 2002; Harris & Zaritsky 2004, 2009; 
Blair et al. 2009; Lawton et al. 2010), and clustering 
analyses have demonstrated the hierarchical formation 
of star clusters and dust clumps (Bonatto & Bica 2010; 
Kim et al. 2010). 

We build upon this work by examining the resolved 
spatial distribution of dust and stellar emission within 
the MCs, using new far-IR and submillimeter (submm) 
photometry from the HERITAGE survey (Meixner et 
al. 2010; Meixner et al., in prep.), near-IR photometry 
from the Spitzer Surveying the Agents of Galaxy Evolu- 



tion (SAGE, Meixner et al. 2006; SAGE-SMC, Gordon 
et al. 2011) surveys, and the Southern Ha Sky Survey 
Atlas (SHASSA, Gaustad et al. 2001). Our goal is to 
explore the energy balance and abundance of dust and 
stars throughout the LMC and SMC, and their implica- 
tions for galaxy evolution. 

This paper is organized as follows. We describe the 
data and their processing in the next section, and in Sec- 
tion^ we describe the dust and stellar properties that we 
analyze in the paper. In Sections HJ [5j and[6j we present 
our main results: we construct andinterpret maps of dust 
luminosity, dust /stellar luminosity and mass, and SFR, 
for the LMC and SMC. We focus on the relation between 
dust/stellar ratios and other properties, especially dust 
heating and star formation. We conclude with a sum- 
mary and discussion of our results. 

2. DATA 

We briefly describe below the data used in this paper. 
We refer the reader to Meixner et al. (2010; Meixner et 
al., in prep.) for details about the observations using the 
Herschel Space Observatory (Pilbratt et al. 2010) and 
data reductions. 

HERITAGE is a uniform survey of the LMC, SMC, 
and Magellanic Bridge with the Spectral and Photomet- 
ric Imaging Receiver (SPIRE) at 250, 350 and 500 fim 
(Griffin et al. 2010) and the Photodetector Array Camera 
and Spectrometer (PACS) at 100 and 160 /im (Poglitsch 
et al. 2010). The SPIRE beam sizes are approximately 
18, 25, and 37" at 250, 350, and 500 ^m, and the flux cal- 
ibration uncertainties are at most 15%. The PACS beam 
sizes are 7 and 11" at 100 and 160 /im, and the calibration 
uncertainties are at most 20%. The data were processed 
using the HIPE 7.0 data reduction software (Ott 2010). 
The images were converted from Jy/pixel (PACS) and 
Jy/beam (SPIRE) to MJy/sr. 

To these we add ancillary data, including imaging from 
the MIPS (Rieke et al. 2004; 24, 70, and 160 fim) and 
IRAC (Fazio et al. 2004; 3.6, 4.5, 5.8, and 8 fj,m) in- 
struments on board Spitzer Space Telescope (Meixner et 
al. 2006; Gordon et al. 2011), and Ha imaging from the 
SHASSA (Gaustad et al. 2001) survey. 

The LMC and SMC multi-wavelength photometry are 
at different resolutions. In order to use these data to- 
gether, it is necessary to convolve them to a common 
resolution (that of SPIRE 500 fim) to generate images 
with a common point spread function (PSF). We use the 
convolution kernels of Aniano et al. (2011), which use 
techniques similar to those of Gordon et al. (2008). We 
then subtract a constant background from the images, 
by fitting Gaussian distributions to the flux densities be- 
yond the galaxies; data within ler of the background are 
not used. We subsequently crop the images and align 
them with each other, yielding image arrays of the same 
dimensions, with 14// pixel size (3.4 and 4.1 pc linear 
size for the LMC and SMC). These are the images that 
we use for computing dust and stellar luminosities and 
masses in the following sections. 

3. DESCRIPTION OF GALAXY PROPERTIES 

In this section, we describe the dust and stellar proper- 
ties of the Magellanic Clouds that we will use throughout 
the paper. Our goal is to analyze the relations between 
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these dust and stellar properties within the MCs. All 
of the properties are inferred from the resolved SEDs, 
except for the SFR, which includes Ha data. 

3.1. Dust Properties 

3.1.1. Dust Luminosity 

The 'dust luminosity' is a useful quantity because it 
can be directly inferred from the IR fluxes, and because 
it can be used as a proxy for the obscured star formation 
as well as the temperature of dust grains (e.g., Dale & 
Helou 2002; Draine & Li 2007). 

Following Skibba et al. (2011; hereafter Sll), we mea- 
sure the dust luminosity Ldust as follows: 

-kdust = 4:irD 2 m X f\ dA, (1) 

where D is the distance to the galaxy, f\ is the flux den- 
sity at wavelength A, and the integration is performed be- 
tween IRAC 5.8/zm < A < SPIRE 500/xm. The resolved 
SED is interpolated between the bands over this range, 
and pixels with noisy or missing far-IR data are not in- 
cluded. Unlike in Sll, we have not subtracted the stellar 
emission from the mid-IR bands before performing the 
integration, as they found that the level of contamina- 
tion was negligible for dwarf and late-type galaxies (less 
than a few per cent); the contamination by starlight is 
negligible throughout the LMC and SMC as well. We 
also do not extrapolate beyond 500/im, as this contribu- 
tion is negligible as well (less than 0.5%). Ldust is usually 
dominated by the fluxes between 70/xm < A < 250/im. 

The luminosities are converted from erg/s/sr to 
erg/s/pc 2 , and then we normalize by the solar luminosity 
Lq. By summing over the pixels of these L^ ust surface 
density maps, we obtain total luminosities of 2.5 x 10 Lq 
for the SMC and 2.3 x 10 s L Q for the LMC. These values 
appear to be consistent with the SED analyses of Israel et 
al. (2010) and Meixner et al. (in prep.), though a direct 
comparison is not possible. 

3.1.2. Dust Temperature 

The dust temperatures and masses are estimated from 
fitting the PACS and SPIRE far-IR and submm (100- 
500 ^m) SED with a single-temperature blackbody (i.e., 
a blackbody of temperature T<i us t) modified by an emis- 
sivity law of A~^: 

f x cx A-^B A (T dust ), (2) 

where B\ is the Planck function. The emissivity expo- 
nent at A < 300 /im is taken to be /? = 1.5. At longer 
wavelengths, the emissivity law is allowed to vary so as 
to fit the surface brightnesses measured in the SPIRE 
350 and 500 ^m bands, and typically varies between 
1 < ft < 2 throughout the galaxies. These values of 
the emissivity index are consistent with previous studies 
(Sakon et al. 2006; Paradis et al. 2011; Planck Collabo- 
ration 2011). This model has been developed to produce 
small fractional residuals and to be consistent with con- 
straints from gas-to-dust ratios (see also Galliano et al. 
2011). We use preliminary dust temperature (and mass) 
maps from Gordon et al. (in prep.), and we refer the 
reader to that paper and Gordon et al. (2010) for details. 



The approach here is similar to that of Engelbracht et al. 

(2010) and Sll. 

We note that the dust temperature Td us t and the emis- 
sivity index /3 are not independent; indeed, an inverse 
correlation between T<j us t and /3 has been observed, which 
some have argued is due to physical properties of the dust 
grains (Dupac et al. 2003; Paradis et al. 2011; Juvela & 
Ysard 2012) or to noise (Shetty et al. 2009; Kelly et al. 
2012). 

The modified-blackbody approach closely reproduces 
the 100-500 /xm PACS and SPIRE photometry, but un- 
derpredicts the MIPS 70 pm photometry. The '70 pm ex- 
cess' (i.e. the excess of the 70 pm. flux densities measured 
by MIPS over that predicted from the single-temperature 
modified-blackbody fitting) is briefly discussed in Ap- 
pendix [B] as arising from stochastically heated small 
grains, or from the fact that the single-temperature as- 
sumption is not realistic. 

3.1.3. Dust Mass 

With T dust (and (3) determined from Eqn. (j2j), we ob- 
tain the dust mass Mdust from f\\ 

A47TP 2 

JWdust - a u it* \ \ 6 ) 

Kabs,A47T-0 A (,J dust ) 

where K a bs,A is the mass absorption coefficient at wave- 
length A. This is modeled as a single component with 
an average temperature. A model including a second, 
cold dust component also yields good fits, though recent 
analyses of far-IR SEDs with Hcrschel data have shown 
that assuming a cold dust component yields dust masses 
that are too large based on metallicity constraints (e.g., 
Bot et al. 2010; De Looze et al. 2012; Aniano et al. 2012; 
Galametz et al. 2012). 

To compute Md us t in Eqn. Q3J) , for spherical grains of 
radius a, K abs ,A = 3 Q a bs(a, A)/i (4 ap), where Q abs {a, A) 
is the absorption efficiency and p is the mass den- 
sity of the dust. For submicron-sized silicate grains 
(p « 3gcm~ 3 ) Q a bs(160/xm) sa 5.5 x 10 -4 , we obtain 
^abs,i60Mm ~ 13.75 cm 2 g _1 (Laor & Draine 1993; Li & 
Draine 2001). The masses are computed at A = 160 pm, 
in order to minimize the dependence on temperature and 
emissivity variations, as well as flux uncertainties (Gor- 
don et al., in prep.). In contrast, Sll used 500 /im as the 
reference wavelength; however, the relative spatial dis- 
tribution is not strongly dependent on this choice. The 
absolute values of the masses may be affected by up to 
0.2 dex by the adopted reference wavelength (see Sll), 
though for low-mass irregular galaxies with significant 
'submm excesses' (e.g., Bot et al. 2010; Galametz et al. 
2011) it is accurate to use 160 pm. 

By summing over the pixels of these dust mass (surface 
density) maps, we obtain total masses of 1.1 x 10 6 M 
for the LMC and 1.1 x 10 5 M Q for the SMC. The LMC 
mass estimate is consistent with the value of 1.2 x 10 6 M Q 
obtained by Bernard et al. (2008) and Galliano et al. 

(2011) , and the SMC mass is approximately consistent 
with 3 x 10 5 M Q obtained by Leroy et al. (2007). The 
spatial resolution here is sufficient to yield unbiased total 
masses (see also Galliano et al. 2011; Aniano et al. 2012); 
based on the results of Galliano et al., the mass values 
may be affected by up to 10%. 
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3.2. Stellar Properties 
3.2.1. Stellar Luminosity 

We will use the 3.6 /an and 4.5 jim photometry as an 
indicator of stellar luminosity and mass in this section 
and the next section. Emission at these wavelengths is 
strongly correlated with the spatial distribution of stel- 
lar light and mass, especially that of old stars, but us- 
ing these wavelengths alone could be contaminated by 
nonstellar contributions (Mentuch et al. 2010; Zibctti 
& Groves 2011; Meidt et al. 2012). For example, hot 
dust and the 3.3 /zm PAH (polycyclic aromatic hydro- 
carbon) emission can cause appreciable contamination, 
contributing ~ 20% of the integrated light at 3.6 /urn in 
star-forming regions (Meidt et al. 2012). The near-IR 
ff-band luminosity would be a superior tracer of stellar 
light and mass (Zibetti et al. 2009), for example, but 
such measurements with resolution and relative sensitiv- 
ity similar to the other data used here do not yet exist. 

Here we will attempt to infer a 'bolometric' stellar lu- 
minosity L*, comparable to the dust luminosity L^ ust 
(Eqn. [lj. For most of the nearby galaxies in Sll, 
was essentially estimated from UV to 5 /iin wavelengths, 
allowing us to constrain the distribution of as a 
function of 1,3.6pm and £3.e/£4.5/«n for galaxies similar 
to the LMC and SMC, with similar morphologies and 
masses. We can also estimate similar global quantities 
for the MCs, using the published SEDs (Israel et al. 
2010; Gordon et al. 2011). The LMC and SMC have 
-^*/-^3.6Mm ~ 80 and 130, respectively, which are slightly 
lower values than the comparable galaxies in Sll (which 
had ~ 100 — 300), but their L 3 6 /L i 5fJir[L ratios are sim- 
ilar. At fixed 3.6 /im luminosity, it is unclear whether 
or to what extent is significantly correlated with 
L3.6/ £4. 5pm, so we choose not to include a correction 
based on this. 

For simplicity, we adopt L. t = /£3.6pm, where / is 
a proportionality constant based on the global stellar 
SEDs of the LMC and SMC (i.e., the factor of 80 or 
130), while acknowledging that this approximation intro- 
duces an uncertainty of ~ 25%, based on the distribution 
of £*/-Zj3.6^m of the comparable nearby galaxies in Sll. 
In addition, from the 3.6 /urn residual (point-source sub- 
tracted) images (Meixner et al. 2006; Gordon et al. 2011, 
we estimate that our values are affected up to 10% 
from contamination by non-stellar sources. The resulting 
L* maps have uncertainties of ~ 0.15 dex. 

3.2.2. Stellar Mass 

To estimate the MCs' stellar mass, we will use the cal- 
ibration of Eskew et al. (2012). This calibration converts 
from the 3.6 and 4.5 /im flux densities (/3.6/im and /4.5pm) 
to stellar mass M*, and it is based on a detailed analysis 
of the LMC, for which they find a scatter of approxi- 
mately 30% in mass around the mean relationship: 



m. = io s - 8B /fiL/ 



-1.85 
' 4.5pm 



D 
O05 



Mr. 



(4) 



where the flux densities /3.6pm and /4.5pm are m Jy and 
the distance D is in Mpc. 

As discussed above, using these wavelengths rather 
than near-IR (e.g., ii-band) fluxes yields a potentially 
biased indicator of stellar mass, and may in some regions 



be affected by hot dust and PAH emission (e.g., Zibetti & 
Groves 2011). In particular, Eskew et al. (2012) note that 
regions that deviate significantly from their calibration 
often have young stellar populations (< 300 Myr old) 
and high 8 /im flux. They also acknowledge that the cal- 
ibration could depend on metallicity, but they were un- 
able to explore this issue because there is little variation 
in metallicity in the LMC (Pagel et al. 1978); however, 
constraining any metallicity dependence of the calibra- 
tion is beyond the scope of this paper. Note that Eskew 
et al. adopt a Salpeter initial mass function (IMF), and 
argue that their inferred stellar mass-to-light ratios mod- 
erately disfavor a bottom- light IMF. The reader should 
bear these caveats in mind. 

The resulting stellar mass (-M*) maps have typical val- 
ues of approximately 120 M Q /pc 2 and 70 M Q /pc 2 in 
the LMC and S MC, respectively. (The dust masses 
in Section 3.1.3 are lower, having typical values of ap- 
proximately O.O7M /pc 2 and O.O2M /pc 2 in the LMC 
and SMC, respectively.) By integrating over the stellar 
mass (surface density) maps, we obtain total masses of 
2.0 x 10 9 Mq for the LMC and 3.1 x 10 8 M Q for the 
SMC, consistent with Harris & Zaritsky (2004, 2009). 

3.2.3. Star Formation Rate 

We now use the combination Ha + 24 /im as a star 
formation rate indicator, in order to account for both 
obscured and unobscured star formation (e.g., Calzetti 
et al. 2007; Kennicutt et al. 2009; Leroy et al. 2012). 

We use Ha data from the SHASSA survey (Gaustad 
et al. 2001). We did not attempt to remove contamina- 
tion by [Nil] emission lines, which could contribute more 
than 20% relative to Ha for some objects (Gaustad et 
al. 2001), but is of order 10% or less for faint dwarf and 
late-type galaxies like the MCs (Kennicutt et al. 2008) 
and is not expected to significantly affect the resulting 
SFRs. 

The SHASSA images are not very deep, and their 
beam size is larger than that of SPIRE 500 fim (~ 76"). 
SHASSA uses a Ha filter centered at 6563A with a 32A 
bandwidth. The LMC and SMC are in fields 013 and 
510, respectively, and for the analysis presented here, 
the SHASSA continuum-subtracted map has been used. 
Since the SHASSA images are at poorer resolution than 
that of SPIRE 500 /im, we fit a Gaussian PSF to the 
continuum images, constructed a convolution kernel, and 
convolved the 24 /im images to this resolution, using the 
code designed by Gordon et al. (2008). A common resolu- 
tion is necessary to combine Ha and 24 /im luminosities 
in the SFR calibration. The SHASSA images are in units 
of deci-Rayleighs, which are converted to our flux density 
units (1R = 2.409 x 10~ 7 erg/s/cm 2 /sr at A = 6563A). 

We compute SFRs of the MCs on a pixel-by-pixel basis, 
using the calibration proposed by Calzetti et al. (2010), 
which is calibrated for normal galaxies and H 11 regions: 



SFR(M Q /yv) = 



CW [I(H«) obs + 0l L(24)] 
if L(24) < 4 x 10 42 erg/s, 

C Hq [£(Ha) obs + a 2 L(24)] 

if 4 x 10 42 < 1,(24) < 5 x 10 43 erg/s, 

(5) 

where the luminosities are in units of erg/s, and Ch q = 
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TABLE 1 

Integrated Properties of the Magellanic Clouds 





log L dust 


log M dust 


log M* 


log SFR 




(log Lq) 


(log M ) 


(log M e ) 


(log M yr" 1 ) 


LMC 


8.4 ± 0.15 


6.0 ±0.1 


9.3 ±0.1 


-0.4 ± 0.2 


SMC 


7.4 ± 0.15 


5.0 ±0.1 


8.5 ±0.1 


-1.6 ±0.2 



The quantities are computed by integrating the re solved 
maps described in this paper. For details, see Sections |3.1.1| 
|3.1.3|[3".2.2| and |3.2.3] about the dust luminosity, dust mass, 
stellar mass, and Sl'K, respectively. 

5.45 x l(T 42 (Af /yr)/(erg/s). The coefficients a x = 
0.020, and = 0.031, are tuned to recover the dust- 
obscured SFR, as opposed to emission due to transiently 
heated dust grains (Calzetti et al. 2007). This calibra- 
tion has a scatter of less than 0.2 dex, and it assumes 
a Kroupa IMF (converting from a Kroupa to Salpeter 
implies SFRs larger by a factor of ~ 1.5). Note that our 
SFR maps may be inaccurate in regions with diffuse 'IR 
cirrus' dust emission, were the dust may be only partially 
heated by young stars (see Kennicutt et al. 2009). 

Our approach differs from the one recently proposed 
by Leroy et al. (2012), which applies detailed dust mod- 
els to the far-IR SED to account for emission due to 
the 'cirrus' radiation field, rather than simply using the 
Ha and 24 /xm luminosities. In addition, their star 
formation surface density is averaged over larger areas 
(~ kpc resolution) than our analysis of the Magellanic 
Clouds. The 24 /jm contribution can result in up to 0.2 
dex uncertainty in the total Ha + 24 /jm SFR; however, 
in both the LMC and SMC, we find that the 24 /jm 
contribution is rather small: in Eqn. [5j 80% and 92% 
of the pixels in the LMC and SMC, respectively, have 
L(24 urn) < 4 x 10 42 erg/s. 

The integrated properties of the LMC and SMC are 
summarized in Tabled] 

4. DISTRIBUTION OF DUST LUMINOSITY AND 

MASS 

We begin by analyzing the distrib ution of dust lumi- 
nosity. As stated in Section |3.1.1| these utilize Her- 
schel data, integrating the luminosity at 5.8 //m < A < 
500 /xm. 

In Figure [T] we show the dust luminosity maps of the 
LMC and SMC. Here and throughout the paper, we ap- 
ply a far-IR signal-to-noise cut, following Gordon et al. 
(in prep.). The dust luminosity is typically in the range 
of 4 - 55 L /pc 2 within the LMC and 3-30 L /pc 2 
within the SMC, with the distributions broadly peak- 
ing at Ldust ~ 7 — 11 and 5 — 8 L /pc 2 , respectively. 
The Ldust distributions of the MCs have "tails" at the 
bright end, as can be seen in the bright regions in the 
figure, and the distributions are not affected by the sig- 
nal/noise threshold used to make the masks, as pixels 
with noisy data and poor S/N have been masked out. 
In Appendix [XJ we compare this dust luminosity to the 
total IR luminosity inferred from shorter wavelengths 
(A < 160 /xm), and show that they have only small dif- 
ferences. 

We now can compare the dust mass to dust luminosity, 
using the Mdust /Ldust ratio, which is shown for the LMC 
and SMC in Figure [2] Md us t /Ldust typically ranges from 
0.001 to 0.005 M Q JL Q in both MCs, with a mean of 



0.0026 Mq/Lq. 

The Mdust /Ldust ratio has an inverted distribution 
with respect to that of Td us t, such that it tends to be 
low in regions where Td us t is high, such as within and 
near the bar of the LMC. 

We show this explicitly in Figure [3[ plotting the L dust 
and Tdust distributions of pixels in the LMC and SMC. 
Approximately 1,000,000 and 110,000 pixels were used 
for the LMC and SMC, respectively. We also plot 
L d ust /Mdust versus Td us t in the right panel of FigurepM It 
clearly shows that the correlation between them is much 
stronger than the Ld U st _ Ld us t correlation: the former has 
Spearman rank correlation coefficients of r s = 0.95 and 
0.92 for the SMC and LMC, respectively, while the latter 
has much weaker values of 0.26 and 0.22pM It is simply 
due to the fact that Ldust oc Afdust x T du J^ (as shown 
by the solid lines in Fig. [3]), while the Ldust _ Ld U st cor- 
relation involves the dust column distribution. For ex- 
ample, there are some regions, such as the southern edge 
of the LMC bar, that have warm temperatures but not 
extremely high luminosities; by normalizing Ld us t with 
Mdust, one effectively removes the dependence on dust 
column. Furthermore, note that the slope of the SMC's 
correlation in the right panel is slightly shallower than 
the LMC's, which implies a slightly lower emissivity in- 
dex & in the SMC. 

Because of the expected spatial variation of dust col- 
umn, one would not expect a one-to-one correspondence 
between Ldust and Td us t- This explains the tail towards 
bright Ldust in the distribution in the left panel of Fig- 
ure [3j This can be understood such that bright regions 
(L d ust > 5OL /pc 2 ) are mostly coincident with dusty 
star-forming regions, because L d„ s t, is an indicator of ob 



scured star formation (see Sec. |3.2.3 ), while the fainter 
regions generally probe the diffuse ISM. 

5. DUST VS STELLAR DISTRIBUTION 

In this section, we show and analyze the spatial dis- 
tribution of the LMC and SMC's dust /stellar luminos- 
ity and mass ratios, using the quantities defined in Sec- 
tion |3j inferred from the resolved SEDs. Galaxy SEDs, 
and dust/stellar ratios in particular, are related to other 
properties indicative of a galaxy's evolution, such as 
mctallicity, stellar mass, starlight intensity, and morphol- 
ogy (e.g., Fontanot et al. 2009; da Cunha et al. 2010). 
The dust/stellar luminosity ratio can be physically in- 
terpreted as the amount of emission being reprocessed 
by dust (both small and large grains) relative to the un- 
obscured emission from (young and old) stars, while the 
dust/stellar mass ratio can be interpreted as the relative 
amounts of stellar mass growth and dust production. As 
shown by Sll, dust/stellar luminosity and mass ratios 
are not equivalent, and have substantial scatter between 
them for dwarf and late-type galaxies. 

5.1. Spatial Distribution of Dust/Stellar Luminosity 

We present maps of the L du *JL* luminosity ratios of 
the LMC and SMC in Figure M The dust /stellar lumi- 
nosity ratios exhibit substantial spatial variation within 
both galaxies. 

19 The Spearman rank correlation coefficient may have a value 
between —1 and 1. A positive (negative) value indicates an 
(anti)correlation, and a value of indicates no correlation. 
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Fig. 1. — LMC (left) and SMC (right) Ldust maps, estimated by summing 5.8-500 fim photometry (see text for details), in units of 
Z/0/pc 2 , as indicated by the color bars below the panels. The right ascension and declination coordinates are in units of degrees; the 
coordinates are omitted from subsequent figures. Green regions are masked, because they have no or noisy far-IR data, or because the 
modified blackbody fits are poor. 
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Fig. 3.— LMC (black solid contours) and SMC (red dot- 
ted contours) distributions of Ldust (left, in Lq/pc 2 units) and 
-kdust/^dust ("ght, Lq/Mq) versus T dust . In the right panel, the 
solid and dashed lines through the contours show idust/^dust 



for /9 = 2 and 1.5, respectively (see text for details). 



The values of Ldust/£* typically range from 0.01-0.2 
and 0.01-0.6 throughout the SMC and LMC, respectively. 
These are comparable to the global values estimated by 
Sll, who obtained L dust /L^ ~ 0.1 — 0.5 for nearby metal- 
poor (Z < OAZq) dwarf galaxies. 

Many regions of the MCs that have bright dust lumi- 
nosities or large dust mass, such as 30 Dor and the SMC 
"bar", not surprisingly also have high L dust /L t ratios. In 
contrast, some regions have particularly low ratios, such 
as the LMC bar (see Figs, fl] and El), which has a relatively 
large stellar flux density (Harris &j Zaritsky 2009). Here, 
dark spots can be seen, and are attributed to regions with 
stars and star clusters (which are resolved in spite of the 
fact that the L* maps are degraded to SPIRE 500 /an 
resolution), lacking significant far-IR emission from dust. 

It is important to keep in mind that stars are discrete 
objects, while dust grains have a more continuous spa- 
tial distribution. In addition, most dust is heated by 
starlight in the diffuse radiation field, while dust near 
compact star-forming regions are heated more intensely 
(e.g., Aniano et al. 2012; Foyle et al. 2012). It is beyond 
the scope of this paper to investigate emission from in- 
dividual stars (see e.g. Harris & Zaritsky 2004, 2009), 
though such an approach is complementary to ours. In 
any case, our dust/stellar ratios are not very sensitive to 
the smoothing scale, and in Section [6] we examine star 
formation activity using smoother Ha luminosity maps. 

One might expect that the dust luminosity is roughly 
proportional to the total dust mass multiplied by the 
stellar luminosity (i.e., L^ ust cx x M,j ust ), based on 
the energy balance of the dust between the absorption 
of starlight (which is proportional to x Md us t) and 
the emission of IR photons (which is essentially L dns t). 
To first order, the spatial variation of Ldust/i* roughly 



reflects the spatial variation of A/dust, as shown by the 
weak correlation in Figure [5](cf. Fig. 4 of Sll, for whole 
galaxies). However, this is complicated by the depen- 
dence on the radiative environment, due to the fact that 
the heating of the dust by starlight is determined by the 
intensity and spectral distribution of the starlight (e.g., 
depending on the dust's proximity to star-forming re- 
gions or whether it is heated by old stars or the diffuse 
ISRF) and by the intrinsic absorption and emission prop- 
erties of the dust (which is determined by its composition 
and size distribution). It is important to keep in mind 
that the stellar emission is due to a variety of stellar pop- 
ulations with different ages, and the dust emission is due 
to dust grains of different compositions and sizes. 

In addition, as we have defined L*, it is more sensitive 
to light due to older rather than younger stellar popula- 
tions. The highest L dust /L^, ratios in Figure E] often occur 
near major Hll region complexes, and may De caused by 
dust directly heated by OB stars and generally hotter 
than elsewhere, but may also reflect a poor recovery of 
the stellar luminosity of the luminous OB star clusters 
exciting these regions (i.e., underestim ating ) by using 
only 3.6 and 4.5 /jm emission (see Sec. 3.2.11. 



5.2. Spatial Distribution of Dust/Stellar Mass 

We now show maps of the mass ratio Md us t/Af* in 
Figure [6] These are qualitatively similar to the maps 
of dust/stellar luminosity (Fig. EL implying that the 
mass and luminosity distributions are similar. Note that 
30 Dor has a M<i us t/M* that is larger than average, but 
its L dust /L^ is particularly large. 

Most regions in the MCs have M dust /-M* < 0.001 
though some regions in the LMC and SMC have ratios 
up to approximately 0.004 and 0.002, respectively. (The 
gradient at the northern edge of the LMC is an artifact, 
due to the different background-subtraction performed 
by Gordon et al., in prep.) We find that the Md us t/M* 
distribution peaks at medians of approximately 3 x 10~ 4 
and 4 x 10 -4 for the LMC and SMC, respectively, while 
L(just/i* peaks at values of 0.14 and 0.07. This is con- 
sistent with Sll, who observed that Ldust/^* is slightly 
more metallicity-dependent than Md us t/-M* in nearby 
late-type and dwarf galaxies. 

The integrated dust /stellar mass ratio for the whole 
LMC and SMC are (M dust ) to t/(M*) tot = 5.4 x 10" 4 and 
3.6 x 10 -4 , respective!}^ We will discuss these ratios in 
the context of star formation and gas in the next section. 

We can compare the LMC and SMC's M dust /M„ to re- 
sults for other galaxies (Dunne et al. 2011; Sll; Cortese 
et al. 2012). Most galaxies have 10" 4 < M dust /M* < 
10 -2 , though the precise values depend on the assump- 
tions involved in both mass estimates. For example, 
Cortese et al. (2012) estimate their stellar masses from 
optical rather than near-IR luminosities and colors, and 
assume a Chabrier IMF, resulting in lower stellar masses 
and hence higher Md us t/-Af*. Nonetheless, the LMC and 
especially the SMC are on the low end of the dust /stellar 
mass ratios in these studies: only three (« l/5th) of the 
dwarf/irregular galaxies in Sll (Holberg II, NGC 5408, 
and NGC 3077) have similarly low ratios, though like 

20 To be clear, note that these are ratios of sums, 
(SMdust)/(SAf«), which are not equivalent to the average ratios, 
(M dust /M*), of the pixels. 
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Fig. 4. — LM C (left) and SMC (right) dust/stellar luminosity (unitless) ratio maps. The dust and stellar luminosities are described in 
Sections I3XT1 and I3T2T1T 




Fig. 5. — LMC (black solid contours) and SMC (red dotted con- 
tours) distributions of L^ ust /L t versus Md ust . 

the MCs, NGC 3077 is tidally interacting with its neigh- 
bors, and a significant amount of dust has been stripped 
(Walter et al. 2011). On average, the SMC's dust/stellar 
ratio is nearly a factor of ten lower than that of low-mass 
dwarfs, while the LMC (which is a one-armed spiral) has 
a « 2x lower ratio compared to similar late-type and 
Magellanic spirals. Dunne et al. (2011), Sll, and Cortese 
et al. (2012) show that dwarf and late-type galaxies that 



are not gas-deficient tend to have higher ratios, as do 
galaxies with higher specific star formation rates; how- 
ever, no strong dependence on mctallicity is observed. 
The MCs are consistent with the dependence on gas con- 
tent: they have slightly low (by a factor of two) dust/gas 
mass ratios compared to galaxies with similar metallici- 
ties (Draine et al. 2007; Herrera-Camus et al. 2012) con- 
sistent with their relatively low dust /stellar masses. 

6. STAR FORMATION AND DUST HEATING 

We present SFR maps of the MCs in Figure [7J They 
provide an empirical description of the spatial distribu- 
tion of star formation in the Magellanic Clouds, related 
to the recent work by Bolatto et al. (2011) and Jameson 
et al. (in prep.), and constraints on the galaxies' star for- 
mation histories (Harris & Zaritsky 2004, 2009). These 
SFR maps are also in approximate agreement with the 
spatial distribution of young stel lar objects (YSOs; Carl- 
son et al. 2012). In Section 6.1 we will show how SFR 
is related to dust heating and emission. 

By integrating over these SFR (surface density) maps, 
and by accounting for the different pixel sizes in or- 
der to conserve the total fluxes, we obtain total SFRs: 
0.38 M Q /yr for the LMC and 0.024 M /yr for the SMC. 
The MCs' very different SFRs are not surprising consid- 
ering their stellar masses, which differ by a factor of five; 
moreover, the LMC appears to lie on the 'star-forming 
sequence' of disk-dominated galaxies (e.g., Brinchmann 
et al. 2004; Schiminovich et al. 2007). Note that, if one 
were to integrate L(Ha) and L(24 fim) and then apply 
the SFR calibration, the upper expression of Eqn. [5] (cal- 
ibrated to normal galaxies) yields SFRs approximately 
30% lower, while the lower expression (calibrated to star- 
forming regions) yields SFRs more similar to the above 



Dust and Stellar Emission of the Magellanic Clouds 



9 





10 



R. A. Skibba, C. W. Engelbracht, et al. 



values. 

These global SFRs are within a factor of two of other 
estimates in the literature, using different methods, such 
as YSO counts, stellar population modeling, and far-IR 
emission. Our LMC SFR is larger than that of other 
estimates (Whitney et al. 2008; Harris & Zaritsky 2009; 
Lawton et al. 2010; Sewilo et al., in prep.), who obtained 
SFR re 0.2-0.25 M /yr. On the other hand, our SMC 
SFR is lower than that of Bolatto et al. (2011) and Wilkc 
et al. (2004), who obtained 0.037 and 0.05 M /yr, re- 
spectively. 

Finally, note that SFR(L IR ) (Kennicutt 1998), yields 
much lower values: 0.006 and 0.08 M Q /yr, for the SMC 
and LMC, consistent with the metallicity dependence 
observed by Dommguez Sanchez et al. (2012). These 
low-metallicity galaxies are dominated by unobscured 
star formation: in 'normal' regions (upper expression in 
Eqn. [5| , unobscured star formation accounts for all of the 
SFR m both galaxies, though in 'H n regions' (lower ex- 
pression, with higher L(24/xm)), obscured star formation 
accounts for 97% and 89% in the LMC and SMC, respec- 
tively. For a more detailed study of the SFR within the 
MCs, we refer the reader to Jameson et al. (in prep.). 

The SFRs can be compared to the star formation his- 
tory maps inferred from optical photometry by Harris 
& Zaritsky (2004, 2009) and atomic and molecular gas 
maps (e.g., Bernard et al. 2008; Bolatto et al. 2011). It 
is clear that the SFR is coincident with known Hn re- 
gions and regions with high gas column densities, and by 
comparing to previous results in this paper, we see that 
in some regions the SFR is coincident with dust emission 
(Fig. [I]) and warm dust (Fig. [2| as well. 

6.1. SFR Correlations 

We now examine these correlations directly, by using 
the SFRs of individual (47"64) pixels. SFR versus the 
dust/stellar luminosity (idust/i*) an d dust temperature 
is shown in Figure [8j Because the images have a large 
number of pixels, we indicate their distributions with 
contours. 

Firstly, it is clear that the SMC's median SFR density 
is a factor of two larger than the LMC's, despite the fact 
that its total SFR is considerably lower. This is partly 
due to a projection effect, because the SMC is more in- 
clined (Groenewegen et al. 2000) than the LMC. It may 
also be due to the different properties of star clusters 
and Hn regions in the SMC and LMC (e.g., Lawton et 
al. 2010; Werchan & Zaritsky 2011), and to a selection 
effect, as low-SFR pixels are easier to detect in the LMC. 

Secondly, a correlation between SFR and dust /stellar 
luminosity or mass has been seen for whole galaxies (da 
Cunha et al. 2010; Sll), and we see a similar, though 
weaker, correlation within the Magellanic Clouds. A 
weak correlation with dust temperature is observed as 
well (consistent with Bernard et al. 2008; Planck Collab- 
oration 2011). The correlation with Tdust has a Spear- 
man rank of r s = 0.45 and 0.36 for the SMC and LMC, 
respectively, and the correlation with Ldust/i* has 0.32 
and 0.33. It is interesting that there is a pronounced 
'tail' in these distributions, with stronger correlations 
at high SFR, L dust /L*, and T dust in the LMC, but it 
has long been known that dust grains may be heated by 
stars, contributing to far-IR emission (e.g., Calzetti et 
al. 1995; Kennicutt 1998; Walter et al. 2007). The sub- 
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Fig. 8. — LMC (black solid contours) and SMC (red dotted con- 
tours) SFR versus T(j ust (left) and Z/dust/-^* (right). The contours 
indicate 10, 20, 33, 50, and 75% of the maximum counts, computed 
from the pixels within the LMC and SMC (Fig. I?}. 

stantial scatter around these correlations demonstrates 
the complex relations between star formation and the 
wavelength dependence of dust heating and emission. 

There are some regions with high Ldust /L* that are not 
rapidly forming stars, especially in the LMC, but they 
generally appear to be in the outskirts of star-forming 
regions. The larger spatial extent of the dust emission 
may simply be due to the fact that the mean free path 
of photons that can heat dust is larger. In addition, the 
southern edge of the LMC's bar, for example, has high 
^dust but less star formation, which is likely an indication 
of a diffuse TR cirrus' component of dust emission. 

6.2. Star Formation Efficiencies 



We measured dust /stellar mass ratios in Section 5.2 



By inverting the ratios, this implies stellar/dust mass 
ratios of approximately 1620 and 3350 for the LMC and 
SMC, respectively. These values can be compared to the 
galaxies' gas/dust ratios, though these have large uncer- 
tainties as well, and there is not yet a consensus in the 
literature on their values. We will take the (atomic plus 
molecular) gas/dust ratios to be (M gas ) to t/(M dust )tot ~ 
340 ± 40 for the LMC (Bernard et al. 2008; Galliano 
et al. 2011; see also Meixner et al. 2010; Roman-Duval 
et al. 2010) and pa 900 ± 150 for the SMC (Bot et al. 
2010b; see also Gordon et al. 2009; Leroy et al. 2011)F] 
with approximate uncertainties based on the range of the 
gas/dust estimates in these papers. These yield star/gas 
mass ratios of 5.5 ± 0.8 and 3.1 ± 0.6 for the LMC and 
SMC, respectively. (In contrast, the MWs star/gas ratio 
is approximately 10.) 
If we were to define a mass-based 'star formation ef- 

21 Note that the gas masses used for these gas/dust ratios include 
a factor of 1.36 to account for the mass of helium (Leroy et al. 2011; 
Galliano et al. 2011). 
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ficiency', SFE m , as the mass ratio of stars/ (stars+gas), 
which is a common definition on the scale of molecu- 
lar clouds (Krumholz & Tan 2007), then we would have 
SFE m!L Mc ~ 0.85 ± 0.15 and SFE m! g M c ~ 0.76 ± 0.20. 

Furthermore, the LMC has a lower gas mass fraction: 
M sas /M* re 18±3%, versus 32±6% for the SMC (though 
not as low as 10% in the MW). This is consistent with the 
fact that more massive galaxies tend to have lower gas 
fractions (Catinella et al. 2010; Masters et al. 2012). An- 
other definition of star formation efficiency of galaxies is 
SFE t = SFR/M gas , whose inverse is the gas consump- 
tion time-scale; this SFE only weakly correlates with 
stellar mass (Bothwell et al. 2009; Huang et al. 2012; 
or with dark matter halo mass, in Dave et al. 2012). By 
this definition, we have SFE t LMC X0 _9 0±01 yr _1 and 
SFE t!S MC ~ 10- 9 - 6±01 yr _1 . SFE t can be interpreted as 
the rate at which available gas is converted into stars, 
while SFE m quantifies the relative amount of gas con- 
verted. Interestingly, the LMC again appears to have a 
relatively large SFE t for its mass, compared to other low- 
redshift late-type and dwarf galaxies (Huang et al. 2012), 
while the SMC's is more typical (see also Bolatto et al. 
2011). In addition, from its star formation history, the 
LMC's SFR has been increasing more rapidly with time 
than the SMC's, on average (Harris & Zaritsky 2009). 

We end by discussing the specific SFR, defined as 
sSFR = SFR/M*. The sSFR is a useful physical quan- 
tity, because by normalizing by the stellar mass, relations 
with star formation activity can appear clearer. In Fig- 
ure |9j we show the distribution of sSFR versus dust tem- 
perature and the dust/stellar luminosity ratio for pixels 
in the LMC and SMC, analogous to Figure IH (and to 
Fig. 7 in Sll, for whole galaxies). We see that the LMC 
has a wider range of sSFRs than the SMC, and that the 
weak correlation between SFR and Td us t disappears when 
accounting for M*. However, the correlation between 
sSFR and Ldust/L* is strong for both galaxies, consistent 
with Sll: regions with larger dust /stellar ratios tend to 
have higher sSFRs. In addition, the sSFR defines a char- 
acteristic time-scale for star formation: Tgp R = AL/M* 
(Brinchmann et al. 2004). This implies that dusty star- 
forming regions, such as 30 Dor and the SMC wing (see 
Fig. Hi), have short star formation time-scales, on the or- 
der of a Gyr. 

7. SUMMARY AND DISCUSSION 
We now summarize our main conclusions: 

• Utilizing the entire dust SED, including Spitzer 
(SAGE) and Herschel (HERITAGE) observa- 
tions, we have spatially resolved dust luminosities 
throughout the LMC and SMC. Some regions with 
bright Ldust coincide with known star-forming re- 
gions. 

• Dust temperature is more strongly correlated with 
^dust/jWdust than Ldust alone. In addition, we de- 
tect a significant 70 /im excess (with respect to 
modified-blackbody fits at longer wavelengths) in 
both galaxies, especially the SMC, indicating the 
presence of stochastically heated dust and/or that 
the dust is heated to a range of temperatures be- 
cause of a range of starlight intensities. 
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Fig. 9. — LMC (black solid contours) and SMC (red dotted con- 
tours) SF R/ M, versus 7d us t (left) and L^ ust /L t (right), analogous 
to Figure [8] 

• Dust/stellar luminosity and mass ratios exhibit 
substantial variation throughout both the LMC 
and SMC. Bright star- forming regions (with con- 
comitant dust emission) have bright Ldust /L*, 
while regions with many stars, especially the LMC 
bar, have low dust/stellar ratios. 

• We probe the spatial distribution of star formation 
in the LMC and SMC, using Ha emission, and find 
that regions with high SFR correlate with those 
with warm dust temperature and bright Ldust/L*. 

• We compare 'global' properties of the LMC and 
SMC to other galaxies (e.g., Skibba et al. 2011). 
The dust/stellar mass ratios, especially of the 
SMC, are lower than that of similar metal-poor 
dwarf and late- type galaxies. The LMC approx- 
imately lies on the stellar mass-SFR relation of 
disk galaxies, but it has a particularly large star 
formation efficiency (quantified by SFR/M gas or 
M*/(M* + M gas )) compared to other galaxies. 

The LMC and SMC appear to be a pair of interact- 
ing and infalling satellite galaxies (perhaps on their first 
passage through the MW halo) , which is a relatively rare 
phenomenon (Liu et al. 2011; Besla et al. 2012). It is an 
interesting question to ask, to what extent these environ- 
mental processes (interaction and infall) are connected 
with the galaxies' evolution, especially their dust pro- 
duction/destruction and star formation histories. It is 
well known that tidal encounters or 'harassment' between 
close pairs of galaxies can modify their morphologies and 
induce bursts of star formation (Moore et al. 1996; Bar- 
ton et al. 2000), but the LMC and SMC allow a close 
examination of the effects of these processes, including 
the indirect effects on dust heating and emission (both 
in star-forming regions and the diffuse ISM). For exam- 
ple, the LMC has essentially become a one-armed spiral 
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galaxy, in which the arm has ongoing star formation and 
substantial dust emission; its off-center bar has limited 
star formation with a relatively warm dust temperature. 

Finally, our results for the global, and especially re- 
solved, galaxy properties of the Magellanic Clouds can 
potentially be used to constrain models, simulations, and 
SED templates. Semi-analytic models of galaxy forma- 
tion make assumptions about the energy balance and 
spatial distribution of stars and dust (e.g., Lacey et al. 
2008; Somerville et al. 2012), while our results constrain 
the strength and scatter in the correlations between dust 
luminosity (or Ld us t /£* ) and dust temperature and SFR 
(Sectionspjandpl). Furthermore, radiative transfer mod- 
els have alxeadybeen compared to the SEDs of nearby 
galaxies (e.g., Misselt et al. 2001; Jonsson et al. 2010; 
Silva et al. 2011) and could benefit from the additional 
constraints on dust properties provided by Herschel (see 
Section [4]), such as the distributions of dust luminosity, 
temperature, and mass. 

7.1. Differences between the Magellanic Clouds 

Though the Magellanic Clouds are often studied to- 
gether, they have many differences. The SMC is more 
dwarf-like than the LMC, which is classified as a barred 
one-armed spiral galaxy (de Vaucouleurs & Freeman 
1972), and has a much lower metallicity than the LMC 
(Russel & Dopita 1992). The MCs also have differ- 
ent gas/dust ratios (Gordon et al. 2003; Leroy et al. 
2011), submillimeter excesses (Bot et al. 2010; Israel et 
al. 2010), 70 jum excesses (Appendix|Bj, dust production 
from evolved stars (Srinivasan et al~2009; Boyer et al. 
2012; Riebel et al. 2012), star formation histories (Harris 
& Zaritsky 2004, 2009), and proper motions (Kallivayalil 
et al. 2006). 

In Section [5l we found that the SMC has a slightly 
lower dust/stellar luminosity and mass ratio than the 
LMC, on average. The SMC's stellar and dust masses 
are five and ten times lower than the LMC's, respectively. 
The SMC also has a much lower SFR and a lower star 
formation efficiency (SFE) than the LMC. 

Many of these different properties of the galaxies could 
be simply explained by the fact that the SMC is consid- 
erably less massive and more metal-poor than the LMC. 
For example, less massive galaxies are expected to have 
lower SFRs than more massive ones, and metal-poor 
galaxies are expected to have higher gas/dust mass ra- 
tios. Nonetheless, the SMC's higher SFR density (Fig. [8) 
and the LMC's relatively high SFE, may be the result of 
these galaxies' unique formation histories, including tidal 
stripping and triggered star formation. 

7.2. Selected Regions in the MCs 

We now discuss a few selected regions: the bar and two 
well-known star-forming regions, 30 Dor and Nil, in the 
LMC, and the SMC wing. (These regions are labeled in 
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Fig. [2}) 

The LMC's bar has a large stellar luminosity and mass, 
indicated by bright near-IR luminosities. It has a large 
fraction of old stars, with ages > 10 Gyr (Harris & Zarit- 
sky 2009). In addition, its dust mass/luminosity ratio is 
relatively low (Fig. [2]) and dust temperature is relatively 
high; considering its limited recent star formation, this 
suggests that the dust is heated by the ISRF. Interest- 
ingly, the LMC's photometric center, which is centered 
on the bar, is offset by ~ 1 kpc from the stellar and Hi 
kinematic center (Cole et al. 2005). Unlike some barred 
spiral galaxies (see the evolving bar fraction in Sheth et 
al. 2008), the LMC's bar has been an integral part of the 
galaxy for most of its history (Harris & Zaritsky 2009) , 
and this asymmetry in the bar appears to have been the 
result of a recent collision with the SMC (Bekki & Chiba 
2007). 

There are large Hn region complexes near the LMC's 
bar, such as 30 Dor and N44 as well as others further 
away, such as Nil in the northwest. These regions 
have substantial stellar emission as well as dust emis- 
sion, though this is not accompanied by a particularly 
large dust mass (i.e., they have low Md us t /-^dust I see 
Fig. f2b . This suggests that these regions have a rela- 
tive underabundance of dust grains, possibly due to dust 
destruction by stellar winds or shocks (Jones et al. 1994; 
Paradis et al. 2009); however, these star-forming regions, 
especially 30 Dor, have enhanced dust/gas ratios (Par- 
adis et al. 2011), indicating that most of their gas has 
been consumed. 

As discussed above, the SMC wing and bridge appear 
to be tidally stripped and their star formation tidally 
triggered (Harris 2007; Gordon et al. 2009, 2011). In ad- 
dition, the wing has a high gas/dust mass ratio, possibly 
due to dust destruction by a harder radiation field and 
shocks during the tidal interaction (Gordon et al. 2009). 
We find that the wing also has high dust/stellar ratios 
(Figs. [4] and [6]), but these are likely due to star forma- 
tion, not to a significant dust abundance, making the 
wing somewhat similar to Hll regions within the MCs. 
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Fig. Al.— LMC (left) and SMC (right) (left) and SMC (right) (left) and SMC (right) (left) and SMC (right) (left) and SMC (right) 
(left) and SMC (right) (left) and SMC (right) (left) and SMC (right) L dust /L Tm ratio. L d ust includes the entire dust SED, while L TIR 
infers it from shorter wavelengths (A < 160 fim). 



APPENDIX 

A. TOTAL INFRARED LUMINOSITY 

We follow Draine & Li (2007), and use a common calibration of the total IR luminosity of dust, inferred from 8, 24, 
70, and 160 /im photometry: 

Ltir = 0.95(^)8.0 + l.lb{vL v ) M + {vL„)to + {vL„)iw (Al) 

where (vL v ) = (L v )c/X. We have also tested the Dale & Helou (2002) formula, which uses only the MIPS bands, and 
have obtained very similar results. Note that with this definition, Ltir is inferred from wavelengths of A < 160 /im, 
in contrast with Ldust (Eqn. IT]), which is summed over the entire dust SED. 

We show the ratio of the dust luminosity Ldust , which includes longer wavelength data, to this Ltir m Figure |A1| 
It is interesting that the Ltir and Ldust maps are similar, but not identical. Their ratio is not uniformly unity; for 
example, it varies in the SE region of the LMC, and is relatively large in the bright star-forming region 30 Dor. This 
is consistent with ongoing work (Galametz et al., in prep.), showing some differences between IR luminosities inferred 
from Spitzer bands (A < 160 /im) and Herschel bands in regions with very high or low surface brightnesses, though in 
general the submm contribution is relatively small in low-metallicity galaxies. The Ldust /Ltir ratio of the pixels in 
the maps has a mean of 1.02± 0.0001 (with standard deviation a = 0.06) in the LMC and 1.07 ± 0.0002 (with a = 0.05) 
in the SMC. The Ltir luminosity inferred from shorter wavelengths is approximately accurate, but estimates based 
on these wavelengths alone could be slightly biased by different grain size distributions or emissivities. By utilizing 
the entire dust SED, Ldust more accurately accounts for the spatial variations of emission from small and large dust 
grains. 



B. 70 MICRON EXCESS 

Using the single-temperature modified blackbody fits to the 100-500 /im dust SEDs, we can compare the 70 /im 
flux predicted from this to the observed MIPS 70 /jm flux. This ratio is an indicator of the spatial distribution of 
stochastically heated small dust grains (Draine & Li 2001; Galliano et al. 2003), and/or the fact that the dust grains of 
different sizes are heated to a range of temperat ures (rather than a single Td us t) by a range of starlight intensities. We 
show maps of it for the LMC and SMC in Figure Bl A ratio (observed/predicted) of greater than unity is evidence for 
a 70 /im 'excess' (Bot et al. 2004; Bernard et al. 2008). We find that most of both MCs have a significant excess, with 
the LMC and SMC hav ing m ean excesses of approximately 1.4 (rms 0.6) and 2.8 (rms 1.0), respectively. The result in 
the left panel of Figure |Blj is not in disagreement with Bernard et al. (2008), who employed a different definition of 
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Fig. Bl. — LMC (left) and SMC (right) 70 /im excess, /70,mips//70 modBBi computed as the ratio between measured MIPS70 fim flux 
and 70 fim flux predicted from modified blackbody model (fitted to 100-500 fim fluxes). 

the 70 [im excess, with respect to a model that includes very small grain emission, and found a strong excess in the 
LMC in regions to the east and south of 30 Doradus. 

Since we observe /7o,MiPs//70,modBB ratios with a significant excess of unity, this implies that a significant fraction 
of small dust grains are stochastically heated by photon absorption. The larger ratios within the SMC is expected 
because the SMC has a steeper UV extinction law, which likely means that the grains in the SMC are generally smaller 
than those in the LMC. 



